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ABSTRACT 


Methods  to  filter  and  to  analyze  measured  data  in  towing  tanks  are  described  and  the 
numerical  results  are  validated  with  the  generated  data.  A  lowpass  filter  is  used  to 
eliminate  the  noise  of  the  measured  data.  The  methods  of  zero  crossing,  sine  function, 
and  spectral  analysis  are  applied  to  compute  the  amplitudes,  periods,  and  phase  angles. 
The  numerical  results  of  die  measured  data  are  compared  with  those  of  the  different 
methods.  When  the  measured  data  is  near  sinusoidal,  the  results  analyzed  with  the  three 
methods  agree  very  well.  The  computer  program,  HARMON  is  written  to  perform  this 
analysis. 


ADMINISTRATIVE  INFORMATION 

The  work  described  in  this  report  was  performed  by  the  Maneuvering  and  Control 
Division  of  the  Hydromechanics  Department,  Carderock  Division,  Naval  Surface 
Warfare  Center.  The  work  was  sponsored  by  the  Office  of  Naval  Research,  Ship 
Structure  and  Systems  S&T  Division  (Code  331)  under  Funding  Document  Number 
1405WX20612.  The  job  order  number  was  05-1-5600-377^0. 


INTRODUCTION 

When  sinusoidal  water  waves  are  generated  by  the  wave  makers  in  towing  facilities, 
the  measured  forces  or  motion  data  of  the  body  being  tested  in  the  facility  are  expected  to 
be  generally  sinusoidal.  Due  to  the  physics  of  real  waves,  in  addition  to  the  noise  in 
measurement  equipment,  the  recorded  data  are  rarely  sinusoidal.  A  numerical  method  to 
analyze  the  recorded  time-history  data  is  necessary  to  find  the  amplitudes  and  phase 
angles  accurately.  The  present  study  will  discuss  the  filtering  of  recorded  data  with  noise 
and  the  analysis  method.  Three  approaches  will  be  discussed:  the  zero-crossing  method, 
the  method  of  sine  function  and  spectral  analysis. 

When  we  analyze  many  data  sets  at  one  time,  we  must  have  a  quick  way  to  analyze 
these  data  and  to  verify  them.  The  main  advantage  of  having  three  methods  is  that,  by 
comparing  the  analyzed  data  with  these  methods,  we  can  distinguish  good  data  from  bad 
data  easily.  When  the  recorded  data  are  nearly  sinusoidal,  the  analyzed  results  from  these 
methods  are  nearly  identical  within  the  limit  of  numerical  accuracy. 

Generated  data  are  used  to  demonstrate  the  validation  of  each  of  these  methods.  The 
numerical  results  are  presented  in  this  report  for  the  analyses  of  the  forces,  moments,  and 
motions  of  a  measured  data  set. 
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DATA  PROCESSING  METHODS 


(A)  FILTERING 

Recorded  time-history  data  usually  contains  noise  due  to  measurement  equipment  or 
vibration  of  the  towing  carriage.  The  noise  in  the  towing  tanks  is  usually  of  high 
frequency  and  a  lowpass  filter  can  be  used  to  remove  it.  The  relationship  of  the  raw  and 
filtered  data  is  as  follows: 

Y(t)  =  \h{r)X(t-T)dT  (1) 


The  function,  h(t)  is  a  filter  since  the  time-history  raw  data  X(t)  is  converted  into  the 
time-history  output  Y(t).  The  Fourier  transform  of  h(t)  becomes  the  transfer  function, 
which  is  a  periodic  function  of  to  ,  equal  to  2rc/A/ ,  as  described  in  Reference  1 . 
Af  is  the  sampling  interval.  The  filter  function  is  expressed  as  the  complex  Fourier 
coefficient  formula  for  the  analysis  of  the  finite  data 


h„  =  —  \H{jo})ejnaiJ  d(0,  -  N  <n<  N 

2/r  J 


jf  f&i 


(2) 


-ttiht 


where  (2N+1)  is  defined  as  the  number  of  samples  in  the  moving  window. 
When  H(jco)  is  real,  Equation  (1)  simplifies  to 


o 


(3) 


In  Equations  (2)  and  (3),  j  denotes  square  root  of  negative  one. 

A  lowpass  non-recursive  digital  filter  with  no  phase  shift  at  any  frequency  is  a  filter 
used  to  eliminate  high-frequency  noise  while  passing  a  low-frequency  signal  without 
changing  its  phase  angle.  The  desired  transfer  function,  Hd(j<o)  is  shown  in  Figure  I  in 
periodic  form.  The  gain  is  equal  to  one  for  frequencies  smaller  than  or  equal  to  the  cutoff 
frequency,  coe  and  is  zero  elsewhere. 

Equation  (3)  now  gives  the  filter  coefficients 

hn  -  ~[  ji-cos(nfijAr)ifry  +  [o  ■  cos(na)At)d(o]  -  sin(nft>,A>)  ^ 
n  i  •  n  ncoAt 

u  c 
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When  n  becomes  zero,  Equation  (4)  is  expressed  as 


o)c  At 

7t 


(5) 


Since  Hd(jw)  in  Figure  1  causes  sharp  discontinuities  in  the  function,  hn,  the  Hanning 
window  is  applied  to  Equation  (4)  as  follows: 

*(0  =  +  c°s(-^4]  (6) 

2  NAl 

The  final  lowpass  filter  is  h„x(t)  and  is  expressed  as  follows: 

I  tir 

=  - sin(ttO)  A/)[l  +  cos(-— — )]  (7) 

2  nx  NAl 

The  Hanning  window  and  lowpass  filtering  functions  with  and  without  the  Hanning 
function  are  plotted  in  Figure  2. 

The  desired  transfer  function  for  a  highpass  non-recursive  digital  filter  with  no  phase 
shift  is  given  in  Figure  lb  and  Equation  (3)  now  gives  the  filter  coefficients 


hn  =  —  [  fo  •  cos(na/At)cl(D  +  [l  ■  cos(no)Al)do>] 

it  i  J 


xS&i 


sin(ff  jr)  sm(«a?eA/) 
nn  nx 


(8) 


Similar  to  the  case  of  the  lowpass  filter,  the  digital  process  for  the  highpass  filter  is 
multiplied  by  the  Hanning  window,  which  is  Equation  (6).  The  details  of  the  filtering 
process  are  given  in  Reference  1 . 

To  verify  the  filtering  process  numerically,  the  following  mathematical  function  is 
used,  which  is  representative  of  a  sinusoidal  signal  with  a  period  of  3.27  seconds,  with 
some  noise  introduced.  This  signal  is  plotted  in  Figure  3a. 

y(t)  -  A{  sin{(W,/ +^, )  +  A2  sin (a)2i  +  <f>2 )  (9) 


where 


4=5.0,  tu,  =1.92, 


A2=  1.0,  co2  =9.42,  <j>2=^~ 
6  4 


(10) 
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The  cutoff  frequency  is  found  as  1.69  times  the  lower  frequency  by  a  trial  and  error 
method.  When  the  cutoff  frequency  is  1.69  times  the  lower  frequency,  the  filtered  data 
are  close  to  the  exact  values.  The  value  of  N  is  calculated  as  follows: 

4;r 

N  =  -^-+1  (11) 

<OcAt 

If  a  time  step  of  A /  =0.01  is  used,  N  becomes  388.  The  Hanning  window  and  lowpass 
filtering  functions  with  and  without  the  Hanning  function  are  plotted  in  Figure  2.  The  role 
of  the  Hanning  window  is  to  reduce  the  magnitude  of  the  filtering  function.  Equation  (4), 
at  the  ends  of  filtering  time  period. 

The  raw  and  filtered  data  are  shown  in  Figure  3.  If  the  filtering  method  is  effective, 
the  filtered  data  will  match  the  first  term  of  Equation  (9).  The  filtered  data  without  the 
Hanning  window  (Fig.  3b)  show  slightly  smaller  amplitudes  than  the  exact  amplitudes, 
which  are  generated  by  the  first  term  of  Equation  (9).  The  plots  of  Figure  3b  arc 
magnified  in  Figure  4.  Figure  4  shows  that  the  values  between  the  upper  and  lower  peaks 
seem  the  same,  but  the  midpoint  of  each  half  cycle  oscillates,  contrary  to  that  of  the 
filtered  data  with  the  Hanning  window  in  Figure  3  c.  The  agreement  between  the 
amplitudes  of  the  filtered  data  with  the  Hanning  window  and  the  exact  data  is  very  good. 

To  filter  the  data  using  a  highpass  filter,  Equation  (8)  and  the  Hanning  window  are 
used,  which  are  shown  in  Figure  5.  The  highpass  signal  could  also  be  obtained  by 
subtracting  the  lowpass  signal  from  the  original  data.  The  filtered  data  with  and  without 
the  Hanning  window  are  shown  in  Figure  6.  When  the  time  is  less  than  3.94  seconds, 
which  is  the  filtering  interval,  the  filtered  data  are  not  usable  since  the  raw  data  arc 
assumed  to  be  zero  for  time  prior  to  zero  seconds.  This  assumption  means  that  the  filtered 
data  for  time  less  than  3.94  seconds  or  for  time  greater  than  the  total  time  minus  3.94 
seconds  are  not  filtered  correctly.  The  exact  data  in  Figure  6  are  the  “noise”  data  from  the 
second  term  of  Equation  (9).  The  amplitudes  of  the  filtered  data  without  the  Hanning 
window  fluctuate  in  time  as  with  the  lowpass  filter.  The  filtered  data  with  the  Hanning 
window  agree  with  the  exact  data  as  shown  in  Figure  6. 


(B)  SPECTRAL  ANALYSIS 

When  the  measured  data  contains  periodic  signals,  the  amplitude  and  period  of  these 
signals  can  be  computed  through  power  spectral  analysis  via  the  Fourier  transform  of  the 
time-history  data.  If  the  measured  data  are  expressed  as  y(t),  the  Fourier  transform  will  be 

QO 

Hay)  =  \y{t)eia>dt  (12) 
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The  inverse  transform  of  Equations  (12)  is 


y(t) 


03) 


If  Ti  is  the  total  time  for  measuring  the  data,  the  power  spectral  density  is  defined  as 
(Reference  2) 


S(a>) 


lim  — 

rr-»*  Tt 


(14) 


For  the  analysis  of  the  measured  data,  the  fast  Fourier  Transform  (FFT)  is  used  instead 
of  Equation  (12)  or  (13)  for  faster  computing  as  detailed  in  Reference  3.  Among  the 
various  methods  of  FFT  listed  in  Reference  3,  the  subroutine  REALFT  is  used  for  the 
present  study.  The  results  of  REALFT  of  y(t)  will  be  complex  values  as  a  function  of  the 
frequency,  w.  The  square  of  the  magnitude  of  the  complex  values  becomes  the  spectra) 
density  and  the  amplitude  of  the  measured  data  will  be  the  square  root  of  two  times  the 
area  of  the  spectral  density  curve,  shown  in  Equation  (14).  The  FFT  of  a  data  set  is 
symmetric.  REALFT  uses  only  the  positive  range,  so  the  results  of  REALFT  will  be  half 
of  the  total  of  Equation  (14). 

The  power  spectral  density  of  a  sinusoid  with  an  amplitude  of  one  (y(t)=sinco0t  with 
ro„=2.45)  is  plotted  in  Figure  7,  which  is  a  delta  function.  Since  the  base  of  the  area  is  two 
times  the  frequency  step,  the  area  under  the  curve  is  two  times  the  total  power  and  the 
square  root  of  this  area  will  be  the  amplitude  of  the  time-history  data.  The  simple  method 
to  compute  the  area  under  the  curve  is  that  it  is  assumed  to  be  a  triangle.  It  is  found  that 
the  area  of  A  ABC  is  1.57.  The  square  root  of  this  value,  which  is  the  amplitude  of  the 
sine  function,  is  1.25.  The  original  amplitude  of  the  sine  function  is  one.  This  indicates 
that  the  assumption  of  triangular  area  under  the  spectral  density  curve  is  not  correct.  In 
order  to  find  the  accurate  amplitude  of  the  time-history  data,  the  section  AB  in  Figure  7c 
is  assumed  as  a  parabola,  which  is  a  convenient  choice,  as 


/(ra)  =  a{&  -  a),  )2  +  / 


where 


a 


(<i>2  -  ft),)2 


(15) 


(16) 
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The  area  under  the  parabola  is  the  integral  of  Equation  (15)  with  the  integral  range 
between  o)t  and  (02  as 

Area  =  +/*)  (17) 

With  the  substitution  of  the  values  in  Figure  7  into  Equation  (17),  the  area  under  the 
parabola  is  0.524,  which  is  half  of  the  total  area.  The  computed  amplitude  is  1 .02. 


(C)  ZERO-CROSSING  METHOD 

The  zero-crossing  method  can  be  used  to  find  the  amplitude  and  period  of  sinusoidal 
or  quasi-sinusoidal  data.  This  method  involves  finding  the  zero  crossing  of  the  data,  and 
computing  the  area  of  the  positive  and  negative  parts.  The  area  can  then  be  used  to 
calculate  the  amplitude  of  the  sine  curve.  The  area  of  the  half  cycle  of  a  sine  curve  with 
an  amplitude  of  a  and  a  frequency  of  to  is  as  follows: 

sin(eot)dt  =  —  =  A  (18) 

a) 

where  a  is  an  amplitude  and  A  is  the  area  under  the  half  cycle  of  a  sine  curve.  The  half 
period  will  be  the  time  interval  between  two  consecutive  times  where  the  magnitudes  arc 
up  or  down  crossing  zero.  By  averaging  the  times  between  successive  zero  crossings,  we 
can  obtain  an  average  half  period  of  the  filtered  data. 

To  find  the  amplitude  of  the  filtered  data,  the  areas  of  the  positive  and  negative 
magnitudes  for  half  cycles  as  shown  in  Figure  8  are  computed  and  the  average  of  the 
absolute  values  of  these  areas  will  be  the  area  of  half  cycle  of  a  sinusoidal  curve.  Then, 
the  amplitude  will  be  computed  by  equating  the  average  area  with  Equation  (18)  as 
follows: 


ft  la i 

Sa 


tw(4+k) 

a  =  v  1  LIE  (19) 

4 

The  present  method  will  be  numerically  validated  with  the  generated  data.  For  the 
single  wave,  the  results  of  the  first  term  of  Equation  (9)  are  used  as  the  generated  data  for 
the  validation.  The  computed  amplitude  is  5  as  given  in  Table  1.  In  Table  I ,  the  results  of 
other  methods,  which  will  be  described  later,  are  also  included. 
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(D)  METHOD  OF  SINE  FUNCTION 
(i)  SINGLE-FREQUENCY  WAVE 

When  the  frequency  of  the  filtered  data  is  known,  a  new  method  to  find  the  amplitudes 
and  phase  angles  of  the  measured  data  will  be  given  here  for  the  single-frequency  data, 
referred  to  as  the  “method  of  sine  function”.  Let  the  measured  data  be  expressed  for  a 
single  wave  as  follows: 


y  =  ^sin((ot  +  <p)  (20) 


where  A  is  the  amplitude,  to  the  frequency,  t  is  time,  and  4>  is  the  phase  angle.  Equation 
(20)  can  be  expanded  as  follows: 


y  =  v4sina#cos^  +  ^costwtsin^  =  Ac  sinev  +  A,  costuf 

(21) 

where 

Ac  =  Acostf;  As=Asixuf>\  A  =  <Jac2+A/ 

(22) 

and 

(a  \ 

<f>  =  tan  1  —  i 

UJ 

(23) 

To  compute  the  amplitudes,  Ac  and  As  in  Equation  (21),  we  substitute  two  measured 
data  points  at  the  times  of  fr  and  ti  into  Equation  (21). 

y^i)  =  Ac  sin  +  As  cos  eof, 

(24) 

y(/2)  =  Ac  sin*; +  As  cos aA2 

(25) 

The  amplitudes,  Ac  and  As  in  Equations  (24)  and  (25)  can  be  computed  easily  by  solving 
these  equations  simultaneously.  However,  when  the  values  of  the  sine  and  cosine 
functions  become  small,  the  numerical  approach  will  cause  an  error  from  dividing  by  a 
small  number.  When  the  argument  of  sine  and  cosine  functions  is  45  degrees,  the  values 
of  the  sine  and  cosine  functions  are  both  maximized.  To  minimize  error,  the  argument  of 
the  sine  and  cosine  functions  in  Equation  (21)  should  be  as  follows: 

cot  =  —  +  m=0,  1,2, ....  (26) 

4 
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To  have  the  values  of  the  sine  and  cosine  functions  in  Equation  (21)  maximal,  the 
measured  data  at  the  following  values  of  time  should  be  used 


ax  1 

m  = - 

n  4 


(27) 


Other  values  for  t  can  also  be  used  as  long  as  the  values  of  sinrot  and  coscot  are  not 
small.  Equation  (27)  is  a  convenient  choice,  but  if  the  precise  value  cannot  be  obtained, 
this  is  not  important.  The  values  of  times  t|  and  t2  in  Equations  (24)  and  (25)  will  be  the 
values  of  t  in  Equation  (27)  for  m  equal  to  zero  and  integers.  In  general,  the  time  in 
Equation  (27)  docs  not  exactly  match  the  time  of  the  recorded  data  because  it  is 
dependent  on  the  time  step  of  that  particular  data  set.  The  time  step  used  should  be  the 
closest  time  step  near  the  computed  time  of  Equation  (27).  There  are  two  solutions  of 
Equation  (21)  for  each  cycle  of  the  measured  data  and  therefore,  there  are  twice  as  many 
solutions  as  there  are  cycles  of  the  measured  data.  The  average  of  these  amplitudes  and 
phase  angles  will  be  the  final  values  for  the  amplitude  and  phase  angle. 

The  present  method  is  numerically  validated  with  the  generated  data.  For  the  single 
wave,  the  results  of  the  first  term  of  Equation  (9)  are  used  as  the  generated  data  for  the 
validation.  The  total  number  of  data  points  is  9000  with  a  time  step  of  0.01  second  and 
with  the  value  of  ft)t,  in  Equation  (9),  m  is  found  to  be  27  different  values  in  this  case 
from  Equation  (27).  There  are  13  sets  of  Equations  (24)  and  (25).  The  averages  of  the 
solutions  of  these  equations  become  the  final  amplitude  and  phase  angle.  The  computed 
amplitude  is  5.0  and  the  computed  phase  angle  is  30.0  degrees  as  given  in  Table  1. 


(ii)  TWO-FREQUENCY  WAVE 

For  the  two-frequency  wave  the  measured  data  can  be  expressed  as  follows: 

y  =  A,  sin(<y,/  +  ft  )  +  A2  sin {m2t  +  <f>2)  +  Ai  sin[(o>,  -m2)t  +  £,]  (28) 

where  the  subscript  1  represents  the  values  of  the  first  frequency  and  the  subscript  2 
represents  the  values  of  the  second  frequency.  Equation  (28)  will  be  expanded  as  follows: 

y  =  A)c  sin  +  Ah  cos  at,/  +  Alc  sin  a>2t  +  A2t  cos  a)2t 
+  AJt.  sin(ru,  ~a/2)t  +  Ais  cos(ryj  - 0)2 )t 


where 
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A\  -  Vv  +  A 


A,  =  Vv+v 

=  “J  Aic  +  /43j 


=  tan' 


-  tan' 


Fa  x 

_k 

\4t  J 


A  > 

^2c  y 


(30) 

(31) 

(32) 

(33) 

(34) 


=  tan' 


u.. 


(35) 


To  Find  the  amplitudes,  A|C,  A|S)  A2C,  Ais,  A3C  and  A3s  in  Equation  (29),  we  substitute 
the  measured  data  at  the  times,  t],  t2,  t3,  U,  ts  and  t$  into  Equation  (29)  as  follows: 


y(t3)  =  Aic  sinew,/,  +  Au  cosew,/,  ■+■  Alc  sinrw2f,  +  A2,  cosew,/, 

+  A3c  sin(ew,  —  ew2)/,  +  4 3j  costal  -  <w2 )/, 

^(/j)  =  Alc  sineu,/2  +  Als  cosew,/2  +  A2c.  sin  ew2/2  +  A2j  cosew2/2 
+  sin(ey,  -ew2)/2  +  AJa  sin(<w,  -ew2)/2 

y(/j)  -  Alc  sin  £W,/3  +  Au  coseu/j  +  A2c  sinew2/3  +  A2s  cos  ew2/3 
+  Aic  sin(ew,  -ew2)/3  +  A3s  sin(£W,  -<w2)/3 

y(t4)  =  Ale  sinew,/4  +  Au  cosew,/4  +  Alc  sinew 2/4  +  A2,  cos<w2/4 
+  A3c  sin(ew,  -  ew2  )/4  +  Ais  sin(ew,  -  o>2  )/4 

y(fs)  =  A ,t.  sineWj/j  +  Au  cosew,/s  +  A2c  sinew2/s  +  Als  cosew2/5 
+  Aic  sin(£w,  - (o2 )/s  +  Ais  sin(ew,  -<w2)/5 

y(t6)  =  A\c  sinew,/6  +  Au  cosew,/s  +  A2c  sinew2/6  +  A1%  cosew2/6 


+  A3e  sin(ew,  -ew2)/fi  +  A3s  cos(<w,  -t»2)/4 


(36) 


(37) 


(38) 


(39) 


(40) 


(41) 
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The  solution  for  A|C,  A[S,  A20  A2S,  A3C  and  A3S  is  similar  to  the  solution  of  Equations  (24) 
and  (25).  To  find  the  time  for  which  the  values  of  the  sine  and  cosine  functions  are 
maximal,  the  arguments  of  the  sine  and  cosine  functions  should  be  as  follows: 


0){t  = 

n 

-  +  /WJT, 

4 

m=0,  1,2, . 

(42) 

co2t  = 

K 

-  +  nn, 

4 

n— (X  1,  2, 

(43) 

To  compute  the  exact  solution  for  the  integers  m  and  n  of  Equations  (42)  and  (43)  is 
difficult.  An  approximate  method  is  applied  to  the  solution  for  these  two  equations.  First 
the  time  for  the  zero  value  of  the  data  is  found,  and  the  half  periods  of  two  frequencies 
are  divided  by  three  different  intervals  by  avoiding  the  small  values  of  sine  and  cosine 
functions.  There  are  twice  as  many  values  for  m  of  Equation  (42)  as  the  number  of 
Equations  (36)  to  (41 ). 

For  the  validation  of  two- frequency  waves,  the  results  of  Equation  (9)  are  used  as  the 
generated  data.  There  are  21  sets  of  the  solutions  of  Equations  (36)  and  (39).  Since 
Equation  (8)  does  not  contain  the  term  for  the  difference  frequencies,  only  four  equations 
are  used  to  validate  the  method.  Therefore,  there  are  the  5  sets  of  Equations  (36)  through 
(39)  and  the  averages  of  the  solutions  of  5  sets  will  be  the  final  amplitudes  and  phase 
angles.  The  computed  amplitudes  are  5.0  and  1.0,  and  the  computed  phase  angles  are  30 
and  45  degrees  as  given  in  Table  2.  In  Table  2,  the  results  of  spectral  analysis  are  also 
included  for  comparison. 

Since  the  measured  data  for  two-frequency  forces  are  not  available,  these  data  are 
generated  with  the  following  equation 

y(t)  =  A,  sinfffljf  +  4 )  +  4  sin(w2t  +  4 ) 

.  (44) 

+  A3  sin[(<u,  -  a)2  )t  +  41  +  4  sin(a>3/) 

where 

4  -  10.0,  o){  =  2.827, 

A2  =5.0,  q}2  =2.199, 

4=1.0,  4=36.0° 

4=1,0,  d>3  =  14.137 

In  Equation  (44)  the  first  two  terms  are  the  contribution  of  the  force  of  each  frequency, 
the  third  term  is  the  difference- frequency  force  and  the  last  term  is  a  noise.  The  computed 
amplitudes  and  phase  angles  are  given  in  Table  3.  The  agreement  between  the  exact  and 
computed  data  is  pretty  good.  The  raw,  filtered  and  recomputed  data  are  plotted  in  Figure 
9.  The  agreement  between  the  filtered  and  recomputed  data  is  very  good. 


fa  =30° 
fa ,  =  45° 


(45) 
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(E)  ANALYSIS  OF  REAL  DATA 


The  results  of  the  measured  forces  of  a  ship  model  undergoing  forced  roll  are 
presented  in  Figures  10  through  13.  Each  figure  consists  of  four  different  plots:  (a)  raw 
data,  (b)  filtered  data,  (c)  filtered  and  predicted  data  from  the  method  of  sine  function, 
and  (d)  spectrum  of  filtered  data.  The  forced  roll  amplitude  is  5  degrees  and  the  period  is 
2  seconds.  The  axial  forces  (Fx)  show  substantial  noise  or  the  presence  of  other 
frequencies  in  the  raw  data,  as  shown  in  Figure  10.  The  computed  amplitudes,  periods, 
and  phase  angles  are  given  in  Table  4.  Since  the  filtered  data  shown  in  Figures  10  and  1 1 
are  sinusoidal,  the  computed  amplitudes  with  the  three  methods  agree  well  as  given  in 
Table  4.  The  areas  of  the  spectra  for  three  forces  are  computed  using  the  parabola  area 
equation  (Equation  (14))  and  the  corresponding  amplitudes  are  the  square  root  of  two 
times  of  the  areas  reported  in  Table  4,  The  axial  and  lateral  force  amplitudes  are 
recomputed  with  the  method  of  sine  function  and  they  are  compared  with  the  filtered  data 
in  Figures  10  and  11,  respectively.  The  agreement  between  the  filtered  and  predicted  data 
is  excellent. 

The  spectrum  of  the  vertical  force  in  Figure  1 2d  shows  two  distinctive  peaks.  This 
indicates  that  the  vertical  force  might  consist  of  forces  at  two  different  frequencies.  The 
spectral  densities  are  maximal  at  the  frequencies  of  3.14  and  6.28  radians  per  second  as 
shown  in  Figure  10.  The  computed  amplitudes  for  these  two  frequencies  are  found  by 
solving  Equations  (36)  to  (39)  without  including  the  difference  terms,  which  are  the  last 
two  terms  in  these  equations.  The  amplitudes  are  1.051  for  the  frequency  of  6.28  and 
0.783  for  the  frequency  of  3.14.  Their  phase  angles  are  -132.94  and  -133.39  degrees, 
respectively.  The  vertical  forces  are  recomputed  with  these  values  and  they  are  compared 
with  the  filtered  data  in  Figure  13.  The  agreement  between  the  filtered  and  recomputed 
data  is  very  good. 

The  results  of  the  measured  moments  of  a  ship  model  undergoing  forced  roll  are 
presented  in  Figures  14  to  16.  The  forced  roll  amplitude  is  5  degrees  and  the  period  is  2 
seconds.  The  lateral  and  vertical  moments  have  some  noise  in  the  signal.  The  computed 
amplitudes,  periods  and  phase  angles  are  given  in  Table  5.  The  recomputed  moments 
with  the  results  of  the  method  of  sine  function  agree  well  with  the  filtered  moments.  The 
spectra  of  all  moments  have  only  one  peak  as  shown  in  Figures  I4d  to  15d. 

The  results  of  the  measured  motions  of  a  ship  model  undergoing  forced  roll  are 
presented  in  Figures  17  to  19.  The  forced  roll  amplitude  is  10  degrees  and  the  period  is  2 
seconds.  The  computed  amplitudes,  periods  and  phase  angles  are  given  in  Table  6.  The 
magnitudes  of  both  heave  and  pitch  motions  are  very  small  compared  with  those  of  roll 
motions.  The  measured  data  of  both  heave  and  pitch  motions  are  not  likely  sinusoidal  and 
there  appears  to  be  an  interference  effect  generated  through  the  interaction  between  heave 
and  roll  motions  or  pitch  and  roll  motions.  Additional  results  of  forced  roll  data  analysis 
are  included  in  Reference  4. 
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CONCLUSIONS 


Methods  to  filter  and  analyze  time-history  data  are  described  and  a  computer  program 
HARMON  is  created  to  find  the  numerical  results.  Based  on  the  present  study  the 
following  conclusions  may  be  drawn: 

1 .  Since  most  noise  of  measured  tank  data  is  high  frequency,  the  lowpass  filter  is  a  good 
choice  to  eliminate  this  noise.  The  highpass  filter  would  be  necessary  to  eliminate  low 
frequency  noise.  The  addition  of  the  Hanning  window  to  the  filter  will  provide  a  better 
match  to  the  original  data. 

2.  When  the  measured  data  are  near  sinusoidal,  the  results  of  the  three  methods:  zero 
crossing,  sine  function,  and  spectral  analysis  in  computing  the  amplitudes,  periods,  and 
phase  angles  agree  well  as  expected.  When  these  results  do  not  agree,  it  would  indicate 
that  the  measured  data  might  not  be  sinusoidal  in  nature. 

3.  Each  method  has  its  own  advantage.  The  method  of  zero  crossing  is  good  for  the 
analysis  of  data  of  single  frequency  for  the  amplitude  and  period.  The  method  of  sine 
function  can  be  applied  to  data  of  known  single  and  two- frequency  data  to  determine 
amplitude  and  phase  angle  for  those  frequencies.  The  spectral  analysis  method  is  good  for 
data  of  single  and  multiple  frequencies,  and  can  be  utilized  even  when  the  frequency  of 
the  signal  is  unknown. 
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(b)  Highpass  transfer  fuction 


Fig.  I .  Lowpass  and  highpass  transfer  functions 
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Fig.  2.  Hanning  window  and  lowpass  filtering  functions 
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Ln 


(a)  Raw  data 


(b)  Filtered  data  without  Hanning  window  and  exact  data 


(c)  Filtered  data  with  Hanning  window  and  exact  data 


Fig.  3.  Raw  data  and  filtered  data  with  lowpass  filter 


Time  (seconds) 


Fig.  4.  Enlarged  plot  of  filtered  data  without  Hanning  window  and  exact  data 
of  Figure  3  b 
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(a)  Hanning  window 


Fig.  5.  Hanning  window  and  highpass  filtering  functions 
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Time  (seconds) 
(a)  Raw  data 


(b)  Filtered  data  without  banning  window  and  exact  data 


Time  (seconds) 

(c)  Filtered  data  with  hanning  window  and  exact  data 


Fig.  6,  Raw  data  and  filtered  data  with  highpass  filter 


S(m).  f(w)  S(w) 


(a)  Sine  wave  of  y=sin{cuot) 


(b)  Spectral  density  of  sino^t 


Fig.  7.  Spectral  density  of  a  sine  function 
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Fig.  8.  Areas  under  half  cycles  of  a  sine  wave  for  zero-crossing  method 
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(d)  Spectrum  of  filtered  data 
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Fig.  9.  Vertical  forces  of  two- frequency  waves  with  frequencies  of  2,827  and 
2A99  radians  per  seconds 
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Fig.  10.  Axial  forces  of  forced  roll  with  the  roll  amplitude  of  5  degrees  and 
the  roll  period  of  2  seconds 
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Fig,  II,  Lateral  forces  of  forced  roll  with  the  roll  amplitude  of  5  degrees  and 
the  roll  period  of  2  seconds 
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(c)  Filtered  and  recomputed  data 


Fig.  12.  Vertical  forces  of  forced  roll  with  the  roll  amplitude  of  5  degrees  and 
the  roll  period  of  2  seconds 
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(c)  Filtered  and  recomputed  data 


(d)  Spectrum  of  filtered  data 


Fig.  13.  Computation  of  the  vertical  forces  with  two  frequencies  of  forced  roll 
with  the  roll  amplitude  of  5  degrees  and  the  roll  period  of  2  seconds 


25 


£ 

s 

x 

2 


-O 

— 

LL 

X 

S 


o 

ED 

«» 

£ 

x 

2 

CO 


20 

15 

10 

-O 

T 

5 

0 

LL 

X 

-5 

2 

-10 

-16 

-.20 

m  vr 

PW  A  l 

^IME  (seconclS) 
(b)  Filtered  data 


?IME  (seconds) 

(c)  Filtered  and  recomputed  data 


Fig.  14,  Axial  moments  of  forced  roll  with  the  roll  amplitude  of  5  degrees  and 
the  roll  period  of  2  seconds 
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(d)  Spectrum  of  filtered  date 


Fig.  15.  Lateral  moments  of  forced  roll  with  the  roll  amplitude  of  5  degrees  and 
the  roll  period  of  2  seconds 
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Fig.  16.  Vertical  moments  of  forced  roll  with  the  roll  amplitude  of  5  degrees  and 
the  roll  period  of  2  seconds 
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SWHE  (ft*sec)  HEVF  (ft)  HEVF  (ft)  HEV  (ft) 


(a)  Raw  data 


TIME  (seconds) 
(b)  Filtered  data 


(c)  Filtered  and  recomputed  data 


Fig.  17.  Heave  motions  of  forced  roll  with  the  roll  amplitude  of  10  degrees  and 
the  roll  period  of  2  seconds 
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{d)  Spectrum  of  filtered  data 


Fig.  18.  Roll  motions  of  forced  roll  with  the  roll  amplitude  of  10  degrees  and 
the  roll  period  of  2  seconds 
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(b)  Filtered  data 


to 

(d)  Spectrum  of  filtered  data 


Fig.  19.  Pitch  motions  of  forced  roll  with  the  roll  amplitude  of  10  degrees  and 
the  roll  period  of  2  seconds 
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Table  I .  Solution  of  the  first  term  of  Equation  (7)  with  three  methods 


Spectral  analysis 

Zero  Crossing 

Sine  function 

Amplitude  (ft) 

5.12 

5.00 

5.00 

Period  (sec) 

3.28 

3.28 

3.28 

Phase  angle  (deg) 

- 

- 

30.00 

Table  2.  Solution  of  Equation  (7)  with  method  of  sine  function  and  spectral  analysis 


Spectral  analysis 

Sine  function 

Amplitude  1  (ft) 

5.12 

5.00 

Amplitude  2  (ft) 

1.02 

1.00 

Period  1  (sec) 

3.28 

3.28 

Period  2  (sec) 

0.67 

0.67 

Phase  angle  1  (deg) 

- 

30.00 

Phase  angle  2  (deg) 

- 

45.00 

Table  3.  Solution  of  Equation  (43)  with  method  of  sine  function  and  spectral  analysis 


Spectral  analysis 

Sine  function 

Amplitude  1 

10.03 

lo.oo" 

Amplitude  2 

5.16 

5.00 

Amplitude  3 

1.03 

1.00 

Period  !  (=T|,  sec) 

2.83 

2.83 

Period  2  (=T2,  sec) 

2.20 

2.20 

Period  3  (=TrT2,  sec) 

0.63 

0.63 

Phase  angle  1 

- 

30,00 

Phase  angle  2 

- 

45.00 

Phase  angle  3 

- 

36.01 
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Table  4.  Computed  amplitudes,  periods  and  phase  angles  of  force  data  for  a  ship 
model  for  forced  roll  with  roll  amplitude  of  5  degrees  and  roll  period  of 
2  seconds 


Spectral  analysis 

Zero  crossing 

Sine  function 

Amplitude  (lb):  Fx 

0.901 

0.890 

0.884 

 Fy  i 

24.985 

24.606 

24.336 

Fz 

1.030 

1.134 

1.118 

Period  (sec):  Fx 

2.047 

1.996 

- 

fy 

2.047 

1.992 

- 

Fz 

2.047 

1.294 

- 

Phase  angle  (deg) :  Fx 

- 

- 

2.32 

fl 

- 

- 

177.07 

Fz 

- 

- 

-153.57 

Table  5.  Computed  amplitudes,  periods  and  phase  angles  of  moment  data  for  a  ship 
model  for  forced  roll  with  roll  amplitude  of  5  degrees  and  roll  period  of 
2  seconds 


Spectral  analysis 

Zero  crossing 

Sine  function 

Amplitude  (ft-lb):  Mx 

15.670 

15.444 

15.322 

My 

1 .587 

1.572 

1.566 

Mz 

18.921 

18.686 

18.671 

Period  (sec):  Mx 

2.047 

1.992 

- 

My 

2.047 

1.994 

* 

Mz 

2.047 

1.994 

* 

Phase  angle  (deg) :  Mx 

* 

- 

177.70 

My 

* 

-178.82 

Mz 

- 

- 

177.34 
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Table  6.  Computed  amplitudes,  periods  and  phase  angles  of  motion  data  for  a  ship 
model  for  forced  roll  with  roll  amplitude  of  10  degrees  and  roll  period  of 
2  seconds 


Spectral  analysis 

Zero  crossing 

Sine  function 

Amplitude:  Heave  (ft) 

0.00355 

0.00385 

0.00389 

Roll  (deg) 

8.16 

8.08 

8.03 

Pitch  (deg) 

0.0920 

0.0929 

0.0907 

Period  (sec):  Heave 

- 

1.13 

- 

Roll 

- 

1.99 

- 

Pitch 

- 

2.08 

- 

Phase  angle  (deg) :  Heave 

- 

- 

21.30 

Roll 

- 

- 

17.10 

Pitch 

- 

- 

147.00 
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